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|. Problematic and choice justification 


Due to climate change and land-use change caused by population and 
economic growth, water resource issues have become more complex. In 
consequence, Improving the reservoir operation for increased efficiency is a 
necessary way of supply management site, which does not require the 
physical development of reservoir. Normally, reservoir operation uses 
upper and lower rule curves to consider the release of water from the 
reservoir responding to downstream demands in long-term operation. 
Optimal rule curves are necessary guidelines in the reservoir operation that 
have been used to assess performance of any reservoir to satisfy water 
supply, irrigation, industrial, hydropower, and environmental conservation 
requirements. The reservoir management agency needs to plan in advance 
the appropriate volume of water in the reservoir (at each time interval) for 
storage and release of water for various purposes, including the 
implementation of the plan as long as the relevant factors in the future 
have not changed from the original. However, if future conditions are 
different from those anticipated in the planning phase, performance may 
differ from planned to minimize water shortage or overflow. The rule curve 
assumptions are based on maintaining the water level (or volume of water) 
in the reservoir to appropriate the changing hydrological situation and 
downstream water allocation. The main purpose is to avoid the risk of 
water shortages and floods in the reservoir and downstream areas. During 
the dry season, reservoirs need to maintain water volume to reduce the risk 
of water levels lower than minimum storage. During the rainy season, the 
reservoir must release water to reduce the storage volume, which can 
support the precipitation and inflow that flows into the reservoir. The 
reservoir rule curves have been improved to provide the optimal solution 
for long-term operation. 


This study proposed a conditional tabu search algorithm (CTSA) to connect 
with the simulation model for searching the optimal reservoir rule curves. A 
minimum average water shortage was used as the objective function for 
the searching procedure. The proposed model has been applied to 
determine the optimal rule curves of the Ubolrat Reservoir in the northeast 
region of Thailand with the historic monthly inflow, future inflow under 
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scenario B2, water demand, hydrologic data, and physical reservoir data. 
Comparison of the conditional genetic algorithm (CGA) and the CTSA was 
shown to demonstrate the effectiveness of the proposed CTSA model. 


ll. Methodology 


1. The Future Inflow 
The development of the optimal future rule curves will use data from 
the future inflow using the PRECIS model (Providing REgional Climates 
for Impacts Studies) considering the effects of climate change. the future 
inflow will be produced using the SWAT hydrological model. PRECIS is a 
regional climate model, based on the development of ECHAM4 model, 
displaying the data as “grid” with high solution of 22 x 22 km?. The data 
recorded during 1997—2014 were used as a baseline to predict those for 
2015 to 2064. 
SWAT (Soil and Water Assessment Tool) is a semidistributed hydrological 
model developed for the measurement of the inflow, sediment, and 
water quality under the climate and land-use changes. SWAT can be 
used to continually measure the daily inflow and define a longer period 
of time in the future. It can also connect and import the spatial data 
from the Geographic Information System (GIS) in order to evaluate the 
inflow. 
The accuracy of the SWAT results can be evaluated by comparing the 
simulated data with that recorded data from the observation station. 
Three variables including R2 (coefficient of determination), RE (relative 
error), and Ens (Nash-Sutcliffe simulation efficiency) were considered as 
the key indicator of the accuracy. In general, SWAT needs to modify the 
value of hydrological parameters for the model calibration and 
validation. In this study, 8 parameter values were used including 
Alpha_BF, Gwqmn, Gw_Revap, Sol_Awc, Epco, Esco, Ch_N2, and 
Gw_delay. Then, the SWAT with adjusted sensitivity parameters was 
optimized. Later, the daily future climate data from PRECIS that had 
been downscaled were classified into 50 years in future periods; the 
processes of model setup are shown in Figure 1. 
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2. Reservoir Operation Model 


Figure 2 
Standard operating rule. 
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The reservoir operation is performed using water usage criteria release, 
operating policies, and reservoir rule curves with monthly data for long- 
term period. A reservoir operation model was constructed on the 
concept of water balance, and it can be used to simulate reservoir 
operation effectively. 

At the point P1, the available water is equal to the water demand, and at 
the point P2, the available water is the total of reservoir capacity and 
water demand. If the available storage is between P1 and P2, water 
release is equal to water demand. If the available water exceeds the sum 
of active storage and water demand, the reservoir will release the excess 
spilled. And: 


D, + Wes = y 


D, for x, < Wyr < Yr +t Da 


for W,, > y, + D, 


r 


\ =y 7 ee 7 A 
D.+W,,—-*, for x,-D, <sW,, <x, 


0 otherwise, 


Was = Syr T Qy = R,, = E, -= DS, (2) 


where R,, is the released water from the reservoir during year v and 
period t(t=1 to 12, representing January to December), D; is the water 
demand of month T; x; is the lower rule curve of month T; y; is the upper 
rule curve of month t, and W,, is the available water calculating by a 
simple water balance, as described in (2). Sy is the stored water at the 
end of month T; Qy, is the monthly inflow to reservoir, E: is the average 
value of evaporation loss, and DS is the minimum reservoir storage 
capacity (the capacity of dead storage). 

In (2) and Figure 2, if available water is in the range of the upper and 
lower rule levels, then demands are satisfied in full. If available water is 
over the top of the upper rule level, then the water is spilled from the 
reservoir in downstream river in order to maintain water level at the 
upper rule level. If available water is under the bottom of the lower rules 
level, a reduction of water release is performed. 
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3. Application of CGA with the Reservoir Simulation Model for Searching 
Rule Curves: 


The CGA and reservoir simulation model for searching the rule curves are 
described in Figure 3 
Figure 3 


Applying CGA and reservoir simulation for searching rule curves. 
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The CGA requires an encoding format to change the decision variables into 
the form of chromosomes. The CGA, which consists of selection, crossover, 
and mutation, is executed. After this stage, the genetic operations will 
create new chromosomes. For this study, each decision variable represents 
the average monthly water storage of the rule curves in the reservoirs, 
which are defined as the upper bound and the lower bound. After the first 
set of chromosomes in the initial population have been calculated (24 
decision variables, which consist of 12 values from the upper bound and 12 
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values from lower bound situations), the released water will be recalculated 
by the reservoir simulation model using these rule curves. Next, the 
released water is used to determine the objective function with the aim of 
assessing the fitness of the GA. After that, the reproduction process will 
create new rule curve values in the next generation. This procedure is 
repeated until the 24 values of rule curves are appropriate. 


the objective function for searching the optimal reservoir rule curves is the 
minimum of the average water shortage (Minjavry)) in million cubic meters 
(MCM) per year, as shown in 


n 


1 
Mina = — ) Sh, (3 
(amy ==), ) 


v=l 


where n is the total number of considered years and Sh, is the water deficit 
during year v (a year that does not meet 100% of the target demand). 


4. Applying Conditional Tabu Search Algorithm for Searching Rule Curves 
The CTSA begins with an initial population {X1, X2, ..., Xn} created 
randomly within the feasible space. With the 24 decision variables (rule 
curve variables for both upper and lower), the feasible solution of the 
iteration ith is represented as X; = [Xi1, X2, ..., Xa]. Then, a set of rule 
curves is used in reservoir simulation, and the released water is 
calculated by the simulation model using these rule curves. Next, the 
released water is used to calculate the fitness function to evaluate the 
feasible solution. The fitness function is the minimum of the average 
water shortage (Z) subject to constraints on the simulation model as 
described in (3). 

Then, the process is continued until the termination criterion is satisfied 
as described in Figure 4. This termination criterion is optimum; it can be 
expressed by a slight change in the fitness values (less than 0.10 MCM). 
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Figure 4 


Applying tabu and reservoir simulation for searching rule curves. 
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5. Illustrative Application 
The Ubolrat Basin is a branch of the Chi Basin located in northeastern 
Thailand (Figure 5). It has an area of about 3,282 km?. The average 
annual rainfall is 1,411 mm, and the mean annual temperature is 27°C. 
The Phong River lies in the middle of the basin. The Ubolrat Dam was 
developed by building an earth core rock fill dam across the Phong River 
with a height of 33 m and crest length of 7,800 m. The normal storage 
capacity and average annual inflow are 2,263 MCM and 2,478.591 MCM, 
respectively. The objectives of the Ubolrat Dam are irrigation, flood 
control, and industrial and domestic water supply. Schematic diagram of 
Ubolrat basin is described in Figure 6. 
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Figure 5 
Location of the Ubolrat Dam. 


101°24°0"E 101°48'0"E 102*12'0"E 102°36'0"E 


Chi Basin 


25 main basins 
of Thailand § 


a — Main stream 
113 — Branch 
A Hydrological station === Basin boundary 


Figure 6 
Schematic diagram of the Ubolrat basin. 
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The study used CTSA in connection with a reservoir operation model to 
find optimal rule curves through the MATLAB toolbox. The optimal rule 
curve can then be applied to an actual scenario depending on whether 
the rule curve can be used to cover every case or event that might occur. 
Thus, the HEC-4 model was used to create the synthetic inflow data into 
the monthly inflows as a synthetic data set of 500 events. This method 
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was based on the actual historic monthly inflow of the Ubolrat Reservoir 
between years 1963 and 2014 (50 years) imported to the HEC-4 model 
to generate the synthetic inflow event. (1 event is a representative 
period of 50 years.) Therefore, the monthly inflow data are 300,000 
values (50 years x 12 months x 500 events). Then, input synthetic inflow 
data were used to assess the efficiency of the new rule curves and 
compare them with the existing rule curves and also between the CTSA 
and CGA model under the same conditions. Further, the new obtained 
rule curves from CTSA and CGA model were used to evaluate with the 
future situation of B2 scenario. The future inflows to reservoir were 
created by SWAT model that considered climate changes. 


Ill. Results and Discussion 


1. Future Inflow into the Ubolrat Reservoir: 
An evaluation on SWAT accuracy used the data found during 1997- 
2014 (18 years; 1997—2008 for calibration and 2009-2014 for 
validation) for Ubolrat Reservoir station. 
The inflow calculated by SWAT and compared with the data from the 
two observed station shows the inflow during the period of model 
calibration and validation; meanwhile, R2, RE, and Ens were 
satisfactory and accurate as the deviation can be accepted. the 
goodness of fit of the data was depicted in Figure 7. 


Figure 7 


The comparison between the runoff from the observed stations and the SWAT result. 
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The inflow at the Ubolrat Reservoir station simulated by SWAT was 
divided into 2 phases: (1) baseline inflow which is the climate and 
spatial data recorded during 1997—2014 and (2) future inflow using 
the climate data from PRECIS model resulted during 2015-2064. The 
inflow analysis indicates that an average volume of the baseline 
inflow was 2,736 MCM and an average volume of the future inflow 
was 4,580.5 MCM. When comparing those two volumes, it was noted 
that the future inflow seems to be increased (1,844.5 MCM or 40.3% 
in 50 future years). Figure 8 illustrates the annual inflow simulated by 
SWAT during 2015-2064, and Figure 9 depicts the comparative result 
between the average baseline inflow and the average 10-year future 
inflow with increased trend in the future period, most of which have 
shown more than 4,000 MCM, except the third period. This result 
shows that the future inflow flows into the Ubolrat Reservoir with 
increased volume. 


Figure 8 


Annual future inflow flowing into Ubolrat Reservoir during 2015—2064. 
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2. Optimal Historic Rule Curves: 
The historic data of inflow, evaporation, water requirement, and 
monthly rainfall were imported for processing in the CGA connected 
to the simulation model and CTSA model, and the optimal rule curves 
were obtained. 


Figure 10 


Optimal historic rule curves of the Ubolrat Reservoir. 
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We observe that the water storage levels of the CTSA and CGA lower 
rule curves are lower than the existing rule curves during the dry 
season (February—June) in order to release more water to reduce 
water scarcity. In the middle of rainy season (August—October), the 
CTSA and CGA upper curves are higher than their existing rule curves 
in order to increase water storage for next dry season. This will help 
alleviate water shortages in the next year. 


. Optimal Future Rule Curves : 
To find the future rule curves, the average monthly inflow for the 
future period 2015—2064 under the B2 scenario was imported into 
the CGA and CTSA model, and then, the optimal future rule curves 
were obtained. Figure 11 shows the optimal upper and lower rule 
curves of the CGA (RC3) and CTSA (RC5) compared with the existing 
rule curves (RC1). 
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Figure 11 


Optimal future rule curves of the Ubolrat Reservoir. 
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4. Performance of Optimal Rule Curves 


Situations of water shortage and excess release of the systems using future inflow. 





Magnitude (MCM/year) Duration (year) 
Situations Rule curves Frequency (times/year) 


Average Maximum Average Maximum 





RC1 (existing) 23.200 660.000 1.167 2.000 
RC2 (CGA) 9.340 412.000 1.200 2.000 
Water shortage RC4 (CTSA) 32.860 418.000 1.417 4.000 
RC3 (CGA) 5.800 290.000 1.000 1.000 
RC5 (CTSA) 32.860 418.000 1.417 





RCI (existing) 2,085.463 4,711.864 24.500 

RC2 (CGA) 2,056.182 4,697.881 24.500 

Excess release RC4 (CTSA) 2,085.965 4,702.331 24.500 
RC3 (CGA) 2,045.096 4,689.417 24.500 

RC5 (CTSA) 2,085.965 4,702.331 24.500 





In the case of future situation (Table), the future rule curves (RC3 and 
RC5) showed the best performance, as indicated by the frequency of 
the water shortage and the average and the maximum magnitudes of 
the water shortages. The future rule curves are more suitable for 





future situations than the existing rule curves and the historic rule 
curve. 


IV. Conclusion: 


1. Is it possible to broaden this subject? 

To search optimal rule curves of the reservoir is a nonlinear 
optimization problem. There are many optimization techniques 
that are applied to connect with the reservoir simulation model 
for searching optimal rule curves such as dynamic programming 
(DP), genetic algorithm (GA), and simulated annealing algorithm 
(SA). Those obtained rule curves are effectively applied for each 
area. However, these new optimization techniques have not been 
applied to find the optimal rule curves like the tabu search 
technique. 

. How about the application in Morocco? 
This method is not hard to use because all we need is historical 
data, spatial data and future expectations from existing studies, to 
start looking for the optimum future curves. It’s something 
necessary and very helpful that help to avoid many catastrophes 
such us floods, as we saw at Casablanca and Tetouan recently. 

. Conclusion 
This study proposed an alternative algorithm for searching optimal 
reservoir rule curves. The conditional tabu search algorithm 
(CTSA) and reservoir simulation model were applied to search the 
optimal rule curves of the Ubolrat Reservoir under historic 
monthly inflow and future inflow under the scenario B2. The 
future inflow and synthetic inflow data of reservoirs were used to 
simulate reservoir system for evaluating situations of water 
shortage and excess release. The results found that the new 
obtained rule curves from CTSA are more suitable for reservoir 
operating than the existing rule curves. The frequency and 
magnitude of water shortage and excess water release for using 
new obtained rule curves are lower than the existing rule curves. 
When comparing the new obtained rule curves from CTSA with 
the rule curves of the CGA method as well as the existing 
simulation method, it was found that these rule curves are similar. 
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The proposed CTSA model is an effective method for application 
to find optimal reservoir rule curves. This reveals that the CTSA 
and GA model with future inflow are effective methods for 
searching optimal reservoir rule curves that are suitable for using 
in the future situations. 





AN 


NS 


ECOLE MOHAMMADIA D'INGÉNIEURS 


V. References 


1. X. Li, S. Gao, P. Liu, and G. Chen, “Dynamic control of flood limited water level for reservoir 
operation by considering inflow uncertainty,” Journal of Hydrology, vol. 391, no. 1-2, pp. 
124-132, 2010. 


P. Liu, L. Li, G. Chen, and D. E. Rheinheimer, “Parameter uncertainty analysis of reservoir 
operating rules based on implicit stochastic optimization,” Journal of Hydrology, vol. 514, 
pp. 102-113, 2014. 


T. Zhao, X. Cai, and D. Yang, “Effect of streamflow forecast uncertainty on real-time 


reservoir operation,” Advances in Water Resources, vol. 34, no. 4, pp. 495-504, 2011. 


E. Fallah-Mehdipour, O. B. Haddad, and M. A. Mariño, “Real-time operation of reservoir 
system by genetic programming,” Water Resources Management, vol. 26, no. 14, pp. 4091- 
4103, 2012. 


S. K. Jain, M. K. Goel, and P. K. Agarwal, “Reservoir operation study of Sabarmati system, 
India,” Journal of Water Resources Planning and Management, vol. 124, no. 1, pp. 31-38, 
1998. 


Electricity Generating Authority of Thailand (EGAT), “Manual for operating the Bhumibol 
and Sirikit Reservoirs,” Bangkok, Thailand, 2004, EGAT Research Report. 


C. Chaleeraktrakoon and A. Kangrang, “Dynamic programming with the principle of 
progressive optimality for searching rule curves,” Canadian Journal of Civil Engineering, 
vol. 34, no. 2, pp. 170-176, 2007. 


F. J. Chang, L. Chen, and L. C. Chang, “Optimizing the reservoir operating rule curves by 
genetic algorithms,” Hydrological Processes, vol. 19, no. 11, pp. 2277-2289, 2005. 


C. H. Yeh, “Applying multi-objective genetic algorithms for planning of flood detention 
dams system,” Taiwan Water Conservancy, vol. 45, no. 2, pp. 70-81, 1997. 


. A. Kangrang, S. Compliew, and R. Hormwichian, “Optimal reservoir rule curves using 
simulated annealing,” Proceedings of the Institution of Civil Engineers-Water Management, 
vol. 164, no. 1, pp. 27-34, 2011. 


. T. Zhang, “A hybrid particle swarm optimization and tabu search algorithm for order 
planning problems of steel factories based on the make-to-stock and make-to-order 
management architecture,” Journal of Industrial and Management Optimization, vol. 7, no. 
1, pp. 31-51, 2011. 





12. 


14. 


15. 


AN 


=Ñ 


ECOLE MOHAMMADIA D'INGÉNIEURS 


M. D. C. Cunha and L. Ribeiro, “Tabu search algorithms for water network 
optimization,” European Journal of Operational Research, vol. 157, no. 3, pp. 746-758, 
2004. 


. F. Glover and M. Laguna, Tabu Search, Kluwer Academic Publishers, Dordrecht, 


Netherlands, 1997. 


U. Faigle and W. Kern, “Some convergence results for probabilistic tabu search,” ORSA 


Journal on Computing, vol. 4, no. 1, pp. 32-37, 1992. 


W. Sa-ngiamvibool, S. Pothiya, and I. Ngamroo, “Multiple tabu search algorithm for 
economic dispatch problem considering valve-point effects,” Electrical Power and Energy 
Systems, vol. 33, no. 4, pp. 846-854, 2011. 


. G. Lacombe, C. T. Hoanh, and V. Smakhtin, “Multi-year variability or unidirectional trends? 


mapping long-term precipitation and temperature changes in continental Southeast Asia 


using PRECIS regional climate model,” Climatic Change, vol. 113, no. 2, pp. 285-299, 2012. 


. Southeast Asia START Regional Center, “Preparation of climate change scenarios for 


climate change impact assessment in Thailand,” Bangkok, Thailand, 2010, Final Report. 


. A. G. Arnold, R. Srinivasan, R. S. Muttiah, and J. R. Williams, “Large area hydrologic 


modeling and assessment pert I: model development,” Journal of American Water 


Resource Association, vol. 34, no. 1, pp. 73-89, 1998. 


. M. T. Vu, V. S. Raghavan, and S. Y. Liong, “Ensemble climate projection for hydro- 


meteorological drought over the river basin in central highland, Vietnam,” KSCE Journal 
of Civil Engineering, vol. 19, no. 2, pp. 427-433, 2015. 


. M. Meaurio, A. Zabaleta, J. Yriarte, R. Srinivasan, and I. Antigtirdad, “Evaluation of SWAT 


models performance to simulate streamflow spatial origin. The case of a small forested 
watershed,” Journal of Hydrology, vol. 525, pp. 326-334, 2015. 


. A. Rittama, “Hedging policy for reservoir system operation: a case study of Mun Bon and 


Lam Chae Reservoirs,” Kasetsart Journal, vol. 43, pp. 833-842, 2009. 


. T. R. Neelakanta and K. Sasireka, “Hydropower reservoir operation using standard 


operating and standard hedging policies,” International Journal of Engineering and 
Technology, vol. 5, no. 2, pp. 1191-1196, 2013. 


. A. Tayebiyan, T. A. M. Ali, A. H. Ghazali, and M. A. Malek, “Optimization of exclusive release 


policies for hydropower reservoir operation by using genetic algorithm,” Water Resources 
Management, vol. 30, no. 3, pp. 1203-1216, 2016. 


18 





ECOLE MOHAMMADIA D'INGENIEURS 


24. A. Kangrang, H. Prasanchum, and R. Hormwichian, “Future runoff under land use and 


climate changes in the Ubolratana Basin, Thailand,” International Journal of Ecology & 
Development, vol. 32, no. 3, pp. 53-66, 2017. 





